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A two- fluid model developed originally to describe wave oscillations in the vacuum arc centrifuge, a cylindrical, 
rapidly rotating, low temperature and confined plasma column, is applied to interpret plasma oscillations in 
a RF generated linear magnetised plasma (WOMBAT) , with similar density and field strength. Compared to 
typical centrifuge plasmas, WOMBAT plasmas have slower normalised rotation frequency, lower temperature 
and lower axial velocity. Despite these differences, the two-fluid model provides a consistent description of 
the WOMBAT plasma configuration and yields qualitative agreement between measured and predicted wave 
oscillation frequencies with axial field strength. In addition, the radial profile of the density perturbation 
predicted by this model is consistent with the data. Parameter scans show that the dispersion curve is 
sensitive to the axial field strength and the electron temperature, and the dependence of oscillation frequency 
with electron temperature matches the experiment. These results consolidate earlier claims that the density 
and floating potential oscillations are a resistive drift mode, driven by the density gradient. To our knowledge, 
this is the first detailed physics model of flowing plasmas in the diffusion region away from the RF source. 
Possible extensions to the model, including temperature non-uniformity and magnetic field oscillations, are 
also discussed. 
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I. INTRODUCTION 

Oscillations have been observed in the electric probe 
signals of many laboratory plasmas, including theta 
pinches fii^ magnetic mirrors, Q machines,^ vacuum 
arc centrifuge^ and helicon plasmas^"— Theta pinches 
and magnetic mirrors behave vastly different from the 
latter because of their high beta and significant mag- 
netic curvature, which result in highly different plasma 
and field geometry. In contrast, Q machine, plasma cen- 
trifuge and helicon plasmas have low beta and negligible 
magnetic curvature, and have similar plasma properties. 
We exploit this similarity by deploying a two-fluid colli- 
sional model, initially developed to describe vacuum arc 
centrifuge plasmas, to describe the configuration and os- 
cillations in a helicon plasma. Our goal is to better de- 
scribe a class of electrostatic oscillations observed, but 
not yet conclusively identified, in the diffusion region of 
helicon plasmas, away from the RF source. 

The helicon plasma we study is provided by the WOM- 
BAT (Waves On Magnetised Beams And Turbulence) 
device, which was designed to provide a plasma environ- 
ment to study wave physics, including chaos, turbulence, 
wave saturation and the related driving mechanisms. 
WOMBAT has been used to study helicon waves fiSiiiii^ 
which have attracted great interest in the past decades 
due to their ability to produce high plasma densitiesJ-S 
Such plasmas also provide a rich environment for the 
experimental study of drift waves and their impact on 
anomalous transport which is important to understand 
in fusion plasmas. Although there have been many pub- 
lications on helicon wave physics, less attention has been 
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devoted to the physics of low frequency oscillations in 
the diffusion region away from the source;^ Light et 
^ji9|20 observed a low frequency electrostatic instability 
only above a critical magnetic field and identified it as 
a mixture of drift waves and a Kelvin-Helmholtz insta- 
bility. Degeling et -'^^ observed relaxation oscillations 
in the kilohertz range that were associated with the var- 
ious types of mode coupling in helicon discharges. More 
recently, a cylindrical linear model which treats global 
eigenmodes?^ has been employed to study drift waves in 
the magnetically confined plasma device VINETA^iiii^ 
however, the model does not involve the E x B rotation 
and neglects ion parallel motion which may influence the 
characteristics of low frequency oscillations;^- Sutherland 
et aP^ observed low frequency ion cyclotron waves that 
were highly localized along the axial center of the WOM- 
BAT plasma device. Spectral measurements revealed a 
four-wave interaction where energy is down-converted to 
the ion cyclotron mode from the helicon pump. 

Low frequency oscillations in helicon plasmas can be 
broadly classified into two types: the Kelvin-Helmholtz 
instability and drift waves. The Kelvin-Helmholtz insta- 
bility is driven by a velocity shear in a mass flow or a ve- 
locity difference across the interface between two fluidsi^ 
This can occur in helicon plasmas, because the streaming 
velocity of the ion fluid is much slower than that of elec- 
trons. The drift wave is a universal instability driven by 
a plasma pressure gradient perpendicular to the ambient 
field. It can arise in fully ionized, magnetically confined 
and low-beta plasmas,-^!, and has been observed in both 
linear and toroidal field geometries i^"— In WOMBAT 
plasmas, the velocity shear is small and the density and 
temperature gradients produce large pressure gradients, 
suggestive of large drive of drift waves. We thus restrict 
attention to this class of modes. 
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In this paper, a two-fluid model, which was developed 
originally for the vacuum arc centrifuge by Hole et a0^^ 
to explain oscillations observed in the density and elec- 
tric potential, is applied to study low frequency oscilla- 
tions observed in WOMBAT. We show that the equilib- 
rium and perturbed density profiles and the space po- 
tential profile from the model and data are consistent. 
The model predicts unstable modes with a similar global 
mode radial structure, and a frequency comparable to 
observed signals. The paper is organised as follows: sec- 
tion II gives a brief description of the experimental setup 
(WOMBAT), section III introduces the two-fiuid model 
and section IV discusses the wave physics revealed by 
this model and the data. Finally, section V discusses the 
possible extensions to this model and presents concluding 
remarks. 



II. EXPERIMENTAL SETUP 




Source region Ditliision chamber 



FIG. 1. A schematic of the WOMBAT (Waves On Magnetised 
Beams And Turbulence) plasma experiment.— 

A schematic}^ of the experiment is shown in Fig. 1. 
The WOMBAT experiment comprises a glass source tube 
and a large stainless steel diffusion chamber, which are 
attached to each other on axis. The source tube is 50 cm 
long and 18 cm in diameter, while the diffusion chamber 
is 200 cm long and 90 cm in inner diameter. A large 
solenoid is employed inside the diffusion chamber to pro- 
vide a steady magnetic field, which can be up to 0.02 T 
and highly uniform along the axis. In this work, the large 
solenoid was used to produce a high density blue core. 
To generate the plasma, a single loop antenna (20 cm in 
diameter, 1 cm wide and 0.3 cm thick) was used in the 
source tube to couple the RF power to the argon gas. The 
7.2 MHz electric current passing through the loop creates 
a time varying magnetic field, which in turn induces az- 
imuthal current in the argon gas and leads to break down 
and formation of the plasma J^i^ The base pressure of the 
diffusion chamber was maintained at 4 x 10^^ Torr by a 
turbomolecular pump and a rotary pump. 

To measure spatial and temporal profiles in the plasma, 
an uncompensated, translating Langmuir probe was in- 
serted radially into the diffusion chamber, 50 cm from the 



source and diffusion chamber interface. The central wire 
of the probe was fed through an alumina support which 
was in turn shielded by a 6 mm diameter grounded steel 
tubing covering the whole extent of the probe length up 
to the probe tips. The removable probe tip was made of 
a 0.2 mm diameter and 8 mm long nickel wire. For the 
present high density plasma discharge in the thin sheath 
regime (the probe sheath is much less than the probe ra- 
dius), the effect brought by RF fluctuations on the float- 
ing potential and the I{V) characteristics was found to 
be negligible)^ Radial translation of the probe was set 
using a computer controlled stepper motor arrangement 
that allowed the probe tip position to be selected with an 
accuracy of a few micrometers. To determine the plasma 
I{V) characteristics, the bias voltage on the Langmuir 
probe was swept between — 50 V and -I- 30 V using 
a Labview program. The plasma density (n^), electron 
temperature (Te), plasma potential (V^) and floating po- 
tential (Vf) were determined from the I{V) characteris- 
tics of the cylindrical probe. 

Typical Langmuir probe measurements of WOMBAT 
argon plasmas are shown in Table I. Other parameters, 
obtained for similar helicon plasmas, '^^'■^^ and measured 
by laser induced fluorescence, include: the ion tempera- 
ture Ti, the bulk rotation frequency ujq, and axial stream- 
ing velocity Vzo ■ The parameter is the argon ion mass 
and the charge number Z is taken to match similar ar- 
gon plasma conditions.'^- The remaining parameters are 
introduced in section III. Finally, plasma parameters for 
the PCEN (Plasma CENtrifuge) device are taken from 
published workj^i^ 



TABLE I. Typical parameters of WOMBAT and PCEN 
(Plasma CENtrifuge). 



parameter 


WOMBAT 


p(-;gj^T33,40 


rii (on axis) 


1.4 X lO^'-* m"^ 


5.2 X 10'^ m 




1.5 eV 


2.9 eV 


Ti 


0.1 eV 


2.9 eV 


mi 


40 amu 


24.31 amu 




0.0185 T 


0.05 T 


Z 


1.0 


1.5 




200 m s"^ 


10* m s-i 


LUO 


1.42 krad s"^ 


184 krad s"^ 


mi 


44.5 krad s"^ 


295 krad s"^ 


= ^ 


0.032 


0.59 




1.72 


1.6 


Bz IE 

_R(characteristic radius) 


0.0235 


0.03 


6 cm 


1.43 cm 



III. PLASMA MODEL 

A. Model assumptions 

The two-fluid model developed by Hole et al^^^^ is 
based on the following plasma assumptions: 
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1. Ions of different charge can be treated as a single 
species with average charge Z . 

2. The plasma is quasi-neutral, giving that rte — Zrii. 

3. The steady-state plasma is azimuthally symmetric and 
has no axial structure. 

4. The effects induced by plasma fluctuations on the ex- 
ternally applied field can be neglected. 

5. Both finite Larmor radius (FLR) and viscosity effects 
are negligible. 

6. For the range of frequencies considered here, the elec- 
tron inertia can be neglected. 

7. The ion and electron temperatures, Ti and Te, are uni- 
form across the plasma column. 

8. The steady-state ion density distribution has a form 
of no = ni{Qi)e~^^/^^ , which is a Gaussian profile. 
Here, ni(0) is the on-axis ion density, and R is the 
characteristic radius at which the density is 1 /e of its 
on-axis value. 

9. The steady-state velocities of ions and electrons 
can be written as Vi = (0, Wjr, v^z) and 
Ve = (0, uje(r)r, Vezir)), respectively, where uji is 
the ion rigid rotor rotation frequency, is the ion 
uniform axial streaming velocity, ujf.{r) is the electron 
rotation frequency, and ^ez(?') is the electron stream- 
ing velocity. While treated in other work*^^, radial dif- 
fusion of both ions and electrons due to electron-ion 
collision is negligible. 

Here, length and time are normalised to R and \/uJic 
respectively, where tOic — ZeBz/rrii is the ion cyclotron 
frequency. A cylindrical coordinate system is developed, 
with (a;, 0, C) = {r/R, 6, z/R) and r = Wict, where 
X and <^ are the normalised radial and axial positions, 
respectively. 



with terms defined as follows: 



Ui 



Vj 



{-Xifii, Xfii, Ui5),Ue 



(x.^e, Xfic, Ue?), 



X ^» / 



■R^ 



h = Ir 



,6 



Bz 7b' 



Here, subscript i and e refer to ion and electron param- 
eters respectively, ip is the normalised radial velocity di- 
vided by X, n is the normalised rotation frequency, m<; 
indicates the normalised axial velocity, A is the ratio be- 
tween ion and electron temperatures, -0 is a convenient 
constant which for A = 1 becomes the square of the nor- 
malised ion thermal velocity, x is a normalised electric 
potential 0, k is the logarithm of the ratio of the ion 
density rii to its on-axis value ni(0), and rig is the ratio 
of the electron-ion collision frequency Vei to its on-axis 
value ;/ei(0). Also, S is the normalised resistivity parallel 
to the magnetic field, where 7]^ is the electrical resistivity 
of a Lorentz gas and is the ratio of the conductivity 
of a charge state Z to that in a Lorentz gas4^ 



C. Steady-state solution 

In cylindrical geometry with a purely axial constant 
field B — (0, 0, Bz), the steady-state solution of this 
model is given by 

xo(^)-xc + [^(l + ^^«) + ^]x^ (5) 



B. Two-fluid equations 

The model comprises the motion and continuity equa- 
tions of ion and electron fluids, written respectively as 



^ + (Ui-V)Ui = -V(ZVx + AVli)+UiX^+&,f (Ue-Ui), 

(1) 



^Z{-VU + Vx) - Ue X <; + (5n,| • (u; - Ue) = 0, (2) 



with 



— = V • Ui + Ui • Vli, 

OT 



dh 
dr 



V • Ue -I- Ue • Vli, 



(3) 



(4) 



Oeo-a,o(l + a,o) + 27/i(A + Z), 



(6) 



where Xc is an arbitrary reference potential. The axial 
current in this model is unconstrained, and can arbitrar- 
ily be set to zero {ui^o — u^^o), consistent with WOM- 
BAT boundary conditions. 

We have compared the steady-state solution to ex- 
perimental data. While there is little spectroscopic 
data available to support the rigid rotor assumption, 
measurements of the fluctuation frequency from probe 
measurements suggest rigid rotation. Figure 2(a) 
shows the equilibrium density proflle of the WOM- 
BAT plasma for two different axial fleld strengths, 
Bz = 0.0185 T and Bz = 0.0034 T, respectively^^ 
Overlaid are best fits to the profile using a Gaussian 
profile (uio = ni{0)e~''^^-'^^ ): these show reasonable 
agreement over the range of the data. The stronger field 
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FIG. 2. Plasma configuration of WOMBAT: (a) Gaussian fit 
of the equilibrium density profile, for Bz = 0.0185 T and 
Bz = 0.0034 T respectively;^^ (b) the radial variation of 
the normalised space potential from the fit procedure (solid 
line), the model (dashed) and the data (dots)j--^ which was 
taken under the conditions of Bz = 0.0068 T, 0.6 mTorr and 
1708 W in the experiment; (c) the radial profile of the electron 
temperature at two field strengths: Bz — 0.0185 T and 
Bz = 0.0034 T.^^ The dashed vertical line in (c) corresponds 
to the radial location of fiuctuation measurements. 



produces better confinement with a smaller characteristic 
radius R. 

Figure 2(b) shows the radial variation of the space po- 
tential in WOMBAT for a slightly different field strength, 
= 0.0068 T.35 Overlaid is a best fit to the 
data using the parabolic potential profile in Eq. (5), 
with the arbitrary reference potential (xc) a-nd gradi- 
ent (G = + riio) + ^) free parameters. Al- 
though the fit, for which G — 0.5, is reasonable, the 
scatter in the data is large. For the same Xc, we have 
over plotted the model potential profile with G = 0.04. 
Within the bulk of the plasma, out to the characteris- 
tic radius r = 0.08 m, model and observed potential 
profiles broadly agree. 

Figure 2(c) shows the radial profile of Tg in the WOM- 
BAT plasmafiS The constant electron temperature of the 
model is an approximation to the radially varying exper- 
imental profile. For the model, we have chosen the Tg 
value at the position at which the frequency was deter- 
mined, r = 7 cm. 



D. Normal mode analysis 

To compute the normal modes of the system, Hole et 
alMi apply a linear perturbation treatment with plasma 
parameters ^ taking the form 



(7) 



Here, s is the perturbation parameter, m the azimuthal 
mode number, kt^ the axial wave number and uj the angu- 
lar frequency. To first order in s the system of Eqs. (1) - 
(4) reduces to 



/ ^{ll,{y) - X[{y)) \ 

^!5^^e-yX[{y) 
\ / 



where A is the matrix 
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Here, uj = uj — mQio — k^u^^ is the frequency in the 
frame of the ion fluid, C = 1 + 2rJio, * = (A + Z)'4), 
y = and we have introduced a new dependent vari- 
able Xi{y), where 

Xi{y) = ^^[ki{y)~xi{y)]. (lo) 

For large axial wavelength modes of the resistive 
plasma column, for which k'^ < S, this model can be 
reduced to a second order differential equation 

(^^^^— ^)L(iV,)[5i(y)]=0, (11) 

where 

^(^c)=.^ + (l-.)^ + (^-^), 

^ {u^-C^){m+^f{y)) ^ mC 
^ U) — milfg + i'^J{y) d) 

and/(y) — i^^e^ with the normalised axial wave number 
F — k^/y/S. For odd TO modes, the boundary conditions 
for Eq. (11) are gi{0) = and gi{Y) = with the 
infinite radius Y representing the edge of plasma column. 
For even to, these conditions become g'i{0) = and 
gi{Y) = 0. We seek unstable solutions to Eq. (11), for 
which > 0. The solutions ui = ± C are stable and 
hence discarded. 



IV. WAVE PHYSICS 

The experimental results indicate that a m = 1 drift 
wave mode dominates the density fluctuations;^ and 
that the mode appears to propagate purely azimuthally 
in direction of the electron diamagnetic drift. There is 
no clear radial component of the mode propagation. The 
observed phase velocity of the mode is approximately 
200 m/s in the region of maximum density perturbation. 



A. Dispersion curve 

As in Hole et alj^^ we have solved Eq. (11) by a shoot- 
ing method. For to = 1 the boundary conditions are 
(71(0) — gi{Y) — 0. As the differential equation 
is homogeneous, the gradient at the edge g'i(Y) is arbi- 
trary: we have chosen g'i{Y) = 1. To solve for given 
_F, we choose a trial Q and march the solution from the 
edge to the core. The complex frequency <D is then ad- 
justed until the on-axis boundary condition is satisfied. 
We commence the procedure at F = 0, for which an 
analytical solution for tD is available. This is found from 
the solution to 



m{uj'^ — C^) mC , , , ^ 

where n is the number of radial nodes in the plasma. 

Figure 3 shows dispersion curve for the plasma con- 
ditions of the PCEN in Table I. We have compared the 
dispersion curve to Hole et al,^^ and found it to be iden- 
tical, thus validating our numerics. For PCEN plasmas, 
the peak of normalised growth rate lj* = 0.39 lies at 
F — 0.3, with normalised frequency uj^ = 0.007: the 
wave is thus near stationary in the frame of the ion fluid. 
The mode crossings located at _F = 0.55 are associ- 
ated with the centrifugal instability. Also shown is the 
dispersion curve for WOMBAT plasma conditions of Ta- 
ble I. For WOMBAT plasmas, the peak growth rate of 
= 0.14 occurs at F = 0.37, for which w'' = 0.44. 
There is no mode crossing in the dispersion curve for 
WOMBAT plasmas. Also, w'' > 0, and so in the labo- 
ratory frame w'' = uj"^ + rnQio + fci;Ui;o, approximately 
uj^ M Lj^ + ilio, such that the frequency of unstable 
waves is always larger than the sum of the plasma ro- 
tation frequency and axial velocity, and the wave prop- 
agates in the -I- 9 direction (electron diamagnetic drift 
direction). We attribute the difference between disper- 
sion curves to the very low ion temperature Ti and slow 
normalised rotation frequency il^o in WOMBAT plasmas, 
compared with PCEN plasmas. 
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FIG. 3. Dispersion curves of PCEN (dashed line) and WOM- 
BAT (solid line), generated by the model based on the con- 
ditions shown in Table I: (a) normalised growth rate a)*; (b) 
normalised frequency (if . (vs normalised axial wavenumber 
F = fc,/V^). 



B. Dispersion curve sensitivity with plasma parameters 

We have examined the sensitivity of the WOMBAT 
plasma dispersion curve with fiio, and Te. For 
riio < 0.05 the dispersion curve changes by less than 
1 %. Figure 4 shows the change in dispersion curve for 
different and T^. Figure 4(a), (b) shows that Gi" and 
a)'" decrease with increasing 3^ : the physical growth rate 
= t^*ajic also obeys this trend while w'' = uJ^oJic 
exhibits the same trend in the range F > 0.5. Also, 
the peak growth rate shifts to lower k^, and so for con- 
stant R the axial wavelength of the most unstable mode 
is increased. In contrast, Fig. 4(c), (d) shows that 
and uj'' increase with increasing T^, and the peak growth 
rate shifts to larger k^. The dependence of growth rate 
with pressure gradient, and frequency with inverse field 
strength is consistent with a resistive drift wavci^ 



C. Wave oscillation frequency with and Te 

In this section we compute the variation in wave fre- 
quency at the maximum growth rate with field strength, 
for two choices of electron temperature: a constant 
Te = 2.13 eV, independent of B^; and a B^ dependent 



Te. As Te data is only available at two field strengths 
(Fig. 2(c)): B^ = 0.0185 T, Te = 3.2 eV; and 
B^ = 0.0034 T, Te = 1.5 eV, we have used lin- 
ear interpolation to compute Te at other field strengths. 
Fluctuation data is available at five field strengths: 
B^ = 0.0129 T, B^ = 0.0146 T, B^ = 0.0162 T, 
B^ = 0.0179 T and B^ = 0.0195 T. We have solved 
the dispersion curve for each case, and selected uj^ corre- 
sponding to the maximum to calculate the frequency 
in the laboratory frame. 

Figure 5 shows the measured and predicted oscillation 
frequencies. For both curves, the predicted and measured 
frequencies differ by a factor of 3.5. There are various 
possible reasons for this gap. The most likely is that the 
model assumes flat Te and T; profiles, whereas WOMBAT 
plasmas exhibit a nonuniform Te profile. It may also be 
the case that the T^ profile is nonuniform. Second, there 
is little data to substantiate the rigid rotor assumption 
of WOMBAT, ■^^ and there may also be slip between the 
rotation of WOMBAT plasmas and the oscillation fre- 
quency. Third, we have neglected the fluctuations in the 
magnetic field, which may affect the dispersion curve. 
Finally, the calculation of Te for each B^ is a linear inter- 
pretation between the two known data, which may bring 
errors and imprecisions into the prediction. 

Figure 5 also reveals that two trends predicted by using 
a constant Te = 2.13 eV and a B^ dependent Te are 
divergent. Predictions for the B^ dependent Te profile 
exhibit the same trend as the data, suggesting the model 
may have correctly captured the Te dependence with Bz ■ 

Finally, to ensure that this model is in principle capa- 
ble to reproduce the data, we adjusted Te to fit the data 
points. We found that a value of Te = 0.5 eV at low 
field dropping to Te = 0.3 eV at high field was able to 
reproduce the observed frequency. 



D. Perturbed density profile 

Figure 6 shows the measured and predicted radial pro- 
files of the perturbed density nn, together with the equi- 
librium density UiQ. The data used here were taken un- 
der the conditions of B^ = 0.01 T, R = 7.6 cm and 
Te ~ 2.46 eV. The linear perturbation treatment gives 
no solution for the absolute magnitude of nn, and so we 
have fitted for the amplitude. Inspection of Fig. 6 sug- 
gests riii has a single peak, and therefore consistent with 
n = of the perturbed mode. Also, the peak of the 
eigenfunction occurs in the region where the equilibrium 
density gradient is large, suggesting that the instability 
has resistive drift- type characteristicsi^l The exact radial 
position of the peak in nn however does not agree. This 
could be due to uncertainty in the core position of the 
plasma, or the assumption of constant Te in the model. 

Finally, to complement our visualisation of 
the mode structure, we have computed the vec- 
tor field of the linear perturbed mass flow by 
TOi(niiUio + nioUii). The perturbed velocity compo- 
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FIG. 4. Results of dispersion curve sensitivity with B^. and Tei (a), (b) scanning results of Bz\ (c), (d) scanning results of Te. 




FIG. 5. Wave oscillation frequency with and Tg 



nentSUa = (x^^ii, xflu, Ui^i),Uel = (x<^ei, xOei, Ue^l) 

and perturbed density riji can be computed from the 
solution of 51(2/) and following equations: 



FIG. 6. Comparison between the measured (dots) and pre- 
dicted (solid) radial profiles of the perturbed density riii, to- 
gether with the equilibrium density riio (dashed). 
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Figure 7 shows the flow vector field at time t = 0. 
As time advances, the mass vector field rotates in the 
clockwise direction, which is the direction of electron dia- 
magnetic drift with into the page. Here, the coordi- 
nate X and y label the cross-section of the plasma, with 
X — rcos9 and y — rsin^. The 9 vector points are 
clockwise with z into the page. The mass flow is zero on 
axis {x = 0, y — 0), consistent with the boundary 
conditions for density. 
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FIG. 7. Vector field plot of the real part of the perturbed 
mass flow for mode (m = l,n = 0) in a cross-section of the 
WOMBAT plasma. 



V. CONCLUSIONS 

In this paper, we employed a two-fluid model, which 
was developed originally for describing the wave oscil- 
lations observed in PCENjHi^i^ to study the low fre- 
quency oscillations in WOMBAT. To ensure that this 



model is consistent with WOMBAT, the measured and 
predicted plasma configurations were compared, includ- 
ing the equilibrium density profile and the space potential 
profile. These show that although the density and space 
potential profiles agree, the temperature profile is not 
flat, as assumed by the model. Next, dispersion curves 
were generated for WOMBAT plasmas. Compared to the 
more rapidly rotating centrifuge plasma, the drift wave 
instability, unstable at larger wavelengths, has a nor- 
malised frequency much larger than the normalised fre- 
quency for the centrifuge. The difference between disper- 
sion curves is principally due to the low rotation speed of 
WOMBAT plasmas compared to those of the centrifuge. 

Our study of the wave oscillation frequency with 
and Te reveals that the measured fluctuation frequency is 
a factor of 3.5 lower than the predicted frequency, and the 
predicted trend of oscillation frequency with B^ for in- 
ferred Tf, matches the data. The discrepancy between the 
measured and predicted fluctuation frequencies may be 
attributable to the limitations of assumed temperature 
uniformity across plasma column. Weaker model limi- 
tations may include the possibility the plasma does not 
rotate rigidly and our neglect of effects induced by plasma 
fluctuations on the externally applied field. Data limita- 
tions include the simultaneous measurements of plasma 
rotation profile and probe frequency, and the measured 
dependency of Tg with B^ ■ 

Finally, we find that the measured and predicted per- 
turbed density profiles have a single peak in the radial 
direction, indicating the perturbed mode has n = 0, and 
the peak position is in the region of maximum equilib- 
rium density gradient. This qualitative agreement con- 
solidates earlier claims that the mode is a drift mode, 
driven by the density gradient of the plasma. To our 
knowledge, this is the first detailed physics model of flow- 
ing plasmas in the diffusion region away from the RF 
source. 
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